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TECHNICAL NOTE 3449 

THEORETICAL ANALYSIS OP INCOMPRESSIBLE FLOW THROIKIH A RADIAL-INLET 
CENTRIFUGAL IMPELLER AT VARIOUS WEIGHT PLOWS 
II - SOLUTION IN LEADING-EDGE REGION BY RELAXATION METHODS 

By James J. Kramer 


SUMMARY 

The detailed solution of the flow aroimd the blade nose of a 48- 
inch-dlameter radial- Inlet centrlftigal Impeller has been obtained by 
relaxation methods for four weight flows. The results are presented In 
a series of figures showing streamlines and relative velocity contours. 
Minimum velocity gradients around the blade nose occurred for the weight 
flow corresponding to a mean angle of attack of -4.6° computed from 
blade speed and an upstream axial-radial velocity for which blade block- 
age has been taken into accomt. A small positive local angle of attack 
seems desirable for blades with rounded leading edges. 


INTRODUCTION 

In reference 1, a method for the solution of the flow through a 
centrifugal impeller was described and applied to an Impeller similar 
to that described in references 2 and 3. The resulting solution was 
then used to check the accuracy of the approximate method described in 
reference 4, especially in the inlet region. More detailed information 
concerning the flow behavior near the leading edge than that obtained 
in the original solution is desirable because a knowledge of the veloc- 
ity gradients in this region is helpful in avoiding boundary- layer sep- 
aration caused by rapidly decelerating flow. The magnitude of losses 
caused by boundary- layer separation at the leading edge is discussed in 
references 5 and 6 for sharp-nosed blades, but the manner in which the 
losses occirr is not discussed. 

The solution was, therefore, refined in the nose region by solving 
for the stream function with the use of a grid of much finer mesh size 
than that used in the original problem. This extension is dependent on 
the original solution because the refinement is obtained in a region 
along whose boundaries the values of the stream f\mctlon are known from 
the original solution. 
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AHALYSIS 

Statement of Problem 

The flow patterns for four weight flows at an opetrating tip speed 
of 700 feet per second for nonviscous, incompressible flow in a radial- 
inlet centrifugal impeller were presented in reference 1. However , the 
information obtained concerning the flow about the leading edge was not 
sufficiently detailed. Therefore, the solutions for the four weight 
flows were refined in the region of the leading edge by solving the par- 
tial differential equation by numerical methods with a grid of consider- 
ably finer mesh. The impeller Investigated, similar to that described in 
references 2 and 3, had a 48-inch diameter and 18 blades. The sharp lead 
ing edge and blunt trailing edge of the test impeller were rounded for 
the theoretical solutions in order to reduce the number of grid points 
necessary to define the leading-edge shape. The blades are shown in the 
radial -tangential plane in figure 1 and the view in the axial -radial 
plane is shown in figure 2. The modified blade coordinates are given in 
table' I of reference 1. The procedure used to refine the flow solution 
around the blade nose is outlined in the following sections . 

Assumptions . - The flow is assumed to be steady, nonviscous, and 
incompressible. The region in which the flow pattern is to be studied 
more carefully is outlined by the bold lines in figure 1. These bound- 
ary lines are chosen sufficiently far from the blade, so that the proc- 
ess of obtaining the refined solution does not significantly affect the 
conditions along the boundaries as defined by the original solution (ref. 
l) . Experience has shomi that the distance between these boundaries and 
the enclosed body should be two to three times the thickness of the en- 
closed body. The assumption is then made that the values of the stream 
functions along these boundaries can be considered known and equal to 
the values obtained in the original solution along these lines . 

Equations . - The definitions, symbols, and equations governing the 
flow are the same as in reference 1, but are repeated herein for conven- 
ience . The solution is obtained on a surface of revolution generated by 
the mean blade-height line. The trace of this surface in the axiatl- 
radial plane is approximated by the following function: 

-0.041456 r 

^ == - 0.40828) (1) 

where r and z are cylindrical coordinates. All symbols are defined 
in the appendix. 

Although the solution is obtained on a surface, a variation in 
stream-sheet thickness (normal to the z-axis) is considered. For this 
problem, the stream-sheet thickness was taken as the blade hei^t (nor- 
mal to the z-axis) and approximated with negligible error by the follow- 
ing equation: 
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b = 0.07208 + 1.01517 


The stream function T 


is defined by the following equations; 


a7 = 

30 = 
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where 9 is the angular cylindrical coordinate of a right-handed sys- 
tem, Wj* and Wg are relative velocity components in the r- and 0- 
directions, respectively, and p is the fluid density. The differen- 
tial equation governing the flow in terms of the stream function under 
the previously stated assumptions is 
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where cc is the angular velocity of the inpeHer and X is equal to 
dr/dz along the surface. This derivative, which can be obtained from 
equatiqn (l), is 


1 _ 0.041456 

^ “ (r - 0.40828)2 

Equation (3) governs the flow and, together with the boundary conditions, 
determines the problem mathematically. 


Boundary conditions . - The boundary conditions are obtained from 
the original solution. The value of the stream function, expressed as 
a dimensionless ratio of the weight flow through a single passage, on 
the blade surface is equal to zeroj the values along the other bound- 
aries are obtained from cross plots of the original data (ref. l) . The 
error involved in reading values from the cross plots was less than 
0.0003 with values of y/M ranging from 0 to 1.0. 


Method of Solution 

Relaxation method . - The solution of equation (3) was obtained by 
means of relaxation techniques (ref. 7). The residuals of the relaxa- 
tion process were reduced to a value indicating unit change in the fifth 
decimal place of the stream function. This degree of accuracy is con- 
sistent with the accuracy of the cross plots used to .obtain the bound- 
ary values . 

GrM. - The grid used in the relaxation process is shown in figure 
3. The density of the grid lines is greatest in the region of the lead- 
ing edge and decreases gradually in both directions. This type of grid 
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assured greatest accuracy in the region of main interest. The ratio of 
the grid spacing in either direction at a point is never greater than 
2.5 and usually less than or equal to 2.0. The intersections of the 
grid lines are called grid points. The grid points in the interior of 
the boundaries are the points at which the numerical solution is ob- 
tained. There were 138 interior grid x>oints for this problem. 

Finite-difference approximation . - The derivatives of equation (3) 
were written in finite-difference foina with a three-point system used, 
that is, with a parabolic variation assumed between grid points. In the 
original solution a five-point system was used. Because of the fineness 
of the grid, a three -point system was considered adequate for this 
problem. 


RESULTS ARD DISCUSSIOH 

The results of the numerical procedure are presented in figures 4 
and 5, which show streamlines (contours of constant 1 /m) and contours 
of constant ratio of relative velocity to tip speed, respectively. Fig- 
ures 4 and 5 are projections on the polar-coordinate plane. The solu- 
tions for four weight flows, the same as those of reference 1, are 
designated as in the following table; 


Case 

Weight flow, 
Ib/sec 

A 

14 

B 

26.25 

C 

32.10 

D 

44 


Streamline^ 

The streamlines, in general, are smooth curves. This result indi- 
cates that a sufficiently fine grid was used in the region of large var- 
iations. 

The stagnation point shifts from the driving face to the trailing 
face as the weight flow increases from 14 to 44 pounds per second. The 
radial position of the stagnation point varies from r ~ 1.05 to 
r •“ 1.028. The relatively close spacing in the radial direction of the 
streamlines for case A is caused by th6 low weight flow. This spacing 
increases as the weight flow increa&es. Although the relative tangen- 
tial velocity just upstream of the blade row is the same for all weight 
flows, the radial spacing of the streamlines changes because of the man- 
ner in which the stream function is made dimensionless in the figure 
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presentation. Ten percent of tlae -weigiit flow through a single passage 
occupies each stream tube for all the examples, so that, in a compari- 
son of the solutions for various weight flows, it should be noted that 
the weight flow through a single stream tube is directly proportional 
to the total weight flow. The streamlines adjacent to the nose curve 
sharply for case A. The c\irvature for case D is not so pronounced be- 
cause of the wider radial spacing of the streamlines. 


Velocities 

Contours of constant ratio of relative velocity to tip speed (tip 
speed of 100 ft/sec for all four weight flows) are shown in figure 5. 

For case A (fig. 5(a)), there is a rapid local acceleration around the 
nose followed by a less rapid deceleration on the trailing face. The 
effect of the eddy attached to the driving face (see ref. l) is mani- 
fested in the almost constant low velocity along the driving face. For 
case D (fig. 5(d)), the very rapid local acceleration followed by the 
less rapid deceleration occurs on the driving face. The peak velocity 
is greater for case D than for case A, and the subsequent deceleration 
is probably more serious than for the lower weight flow because the low- 
momentum air caused by the deceleration aggravates the secondary-flow 
conditions. These secondary flows transport the low-momentum fluid on 
the driving face to the trailing face. This type of motion is discussed 
in more detail in reference 8. 

For case B a slight deceleration occurs on the driving face. As 
pointed out in reference 1, the weight flow of 26.25 pounds per second 
is that for which a mean angle of attack of -4.6° occurs when cOTiqputed 
from the weight flow, rotational speed, and inlet annular area with 
blade blockage taken into account. The mean angle of attack is the 
angle between the tangent to the blade mean line at the leading edge 
and the mean flow direction. The sign convention for angle of attack, 
shown in figure 3, is chosen so that flow directed at the driving face 
results in a positive angle of attack. The weight flow of 26.25 pounds 
per second is also that corresponding experimentally to maximum effi- 
ciency. However, since the leading edge was rounded, the experimental 
results cannot be explained in terms of the exact solution. For the 
solution of case B, the local angle of attack, that is, the angle be- 
tween the tangents to the blade mean line and the stagnation streamline, 
is approximately +6°, whereas for case C the local angle of attack is 
-7°. Thus, the weight flow for zero local angle of attack is about 29 
pounds per second. However, zero local angle of attack is not necessar- 
ily desirable for rounded leading edges such as the one considered here- 
in. In view of the velocity contours of the 26.25- and 32.10-pound- 
per-second cases, a slight positive local angle of attack seems desir- 
able. For a weight flow of 29 pounds per second a greater deceleration 
would probably occur on the driving face than that which occurs for 
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26.25 poutLds per second. This appears to be the case since the leading 
edge is shaped so that at zero local angle of attack the flow on the 
blade surfaces would be roughly symmetrical about the blade mean line. 
Consequently^ it seems advisable to design the leading-edge regions in 
the manner suggested in reference 9 and further discussed in reference 
10. These leading-edge regions are characterized by very little curva- 
ture of the driving face so that flow alined with the driving face would 
produce little or no local deceleration. 


SUMMARY OF RESULTS 

The detailed solution of the flow around the blade nose of a 48- 
inch-diameter radial-inlet centrifugal Impeller was obtained by means 
of relaxation techniques. The results are presented in a series of 
figures showing streamlines and contours of constant resultant velocity 
ratio. The following results were noted: 

1. The stagnation point shifted from the driving to the trailing 
face of the blade as the weight flow increased. The total shift in the 
radial direction was approximately 0.02 foot. 

2. Minimum velocity gradients aromd the blade nose occurred for 
the weight flow corresponding to an angle of attack of -4,6° computed 
from blade speed and an upstream radial-axial velocity for which blade 
blockage has been taken into account. 

3. A small positive loceil angle of attack seems desirable for 
rounded leading-edge blades. 


National Advisory Committee for Aeronautics 
Lewis Flight Propulsion Laboratory 
Cleveland, Ohio, March 4, 1955 
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APPENDIX - SYMBOLS 

The following symbols are used in tliiB report: 
b stream- sheet thickness in z -direction, ft 
M weight flow through single passage, Ib/sec 

r radial distance, ft 

W ratio of relative velocity to tip speed 
w relative velocity, ft/sec 

z axial distance, ft 

6 angular coordinate, radians 

X slope of trace of stream surface in axial- radial plane 
p fluid density, Ib/cu ft 

Y stream ftinction, eq. (2) 

0 ) angular velocity of iurpeller, radlan/sec 
Subscripts: 

r component in radial direction 
t tip 

0 component in tangential direction 
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Hgnre 4. - Continued. StaToenllnes for refined solution In nose region. 
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Ilguro 5. - Ccmtoura of oonrtamt Telocity ratio W for refined •olntlon. 
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(c) Vel^t flow, 32.10 pounds por second (case C) . 


Jlgure 5. - Continued. Contours of constant Telocity ratio V for 

ref ined solution In nose region. } 
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